Adiabatic Preparation of a Heisenberg Antiferromagnet Using an Optical Superlattice 
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We analyze the possibility to prepare a Heisenberg antiferromagnet with cold fermions in optical 
lattices, starting from a band insulator and adiabatically changing the lattice potential. The nu- 
merical simulation of the dynamics in ID allows us to identify the conditions for success, and to 
study the influence that the presence of holes in the initial state may have on the protocol. We also 
extend our results to two-dimensional systems. 



Ultracold atoms trapped in optical lattices offer a 
unique possibility to experimentally explore strongly cor- 
related states of quantum matter. Currently, one of the 
main experimental challenges in this field is the prepa- 
ration of a Heisenberg antiferromagnet (AFM), which 
represents the necessary next experimental step towards 
a true quantum simulator of the fermionic Hubbard 
model pj. 

Although the creation of a fermionic Mott insulator 
(MI) has recently been reported [l, the realization 
of antiferromagnetic order requires temperature and en- 
tropy significantly lower than presently achieved 0, 0] , 
despite many existing proposals for direct cooling within 
the lattice [6| . An alternative to the direct generation is 
to use an adiabatic protocol 0, @1- In such a scheme, it 
is desirable to tune interactions initially to give a ground 
state with very low entropy. Then, they are changed 
slowly, until the Heisenberg Hamiltonian is realized at 
the end. If the process is adiabatic, the entropy will 
stay low and the final state will be the desired AFM. 
The following questions immediately arise: What are the 
conditions to achieve adiabaticity? What occurs if these 
conditions are too restrictive and cannot be met? And, 
how will the protocol be affected by a finite temperature 
and the presence of a harmonic trap? 

In this Letter, we propose a specific adiabatic scheme 
and analyze these issues. Our adiabatic protocol is the 
first to attain an AFM with ultracold fermions within 
feasible timescalcs, even in the presence of experimental 
imperfections. Additionally, we show that it is possible to 
realize antiferromagnetic order on a part of the sample 
in a shorter time than required for the whole system. 
Finally, we simulate the dynamics of holes to demonstrate 
their destructive effect on the AFM and devise a strategy 
to control them. 

The initial ground state of our protocol is a band in- 
sulator (BI), which is transformed first to an array of 
decoupled singlets and finally to the AFM, by adiabat- 
ically changing the depth of two superimposed optical 
lattices. A BI is easier to prepare with low entropy than 
a MI for two reasons. On the one hand, its energy gap 
is given by the band gap, which is much larger than the 



interaction energy (MI gap) and favors a redistribution 
of the entropy towards the surrounding metallic shell [§] . 
On the other hand, the preparation can be done using 
weakly or noninteracting atoms, thereby avoiding the 
long timescales associated with mass and entropy trans- 
port at higher interactions [lo| . 

For the one-dimensional case, we simulate the 
fermionic t — J model with matrix product states 
(MPS) [111. We first identify the adiabatic conditions 
that allow the preparation of the antiferromagnetic state 
in an ideal case with no defects. Second, we study how 
these conditions arc relaxed if one imposes that antifer- 
romagnetic order is only obtained on a subset of fermions 
around the center of the sample. We observe that, when 
restricted to a middle sublattice, adiabaticity is deter- 
mined by an effective gap related to this sublattice and 
not by the gap of the total system. Third, we include the 
presence of holes in the initial state, expected to occur 
in real experiments due to the finite temperature. The 
large initial energy of the holes can in principle destroy 
the AFM as they delocalize inside the sample. We find 
that, if the holes are initially located at the outer part of 
the sample, as expected in an experiment, a tradeoff can 
be reached between the degree of adiabaticity of the pro- 
cess and the distance the holes travel inside the chain, 
so that the antiferromagnetic order is still produced in 
the center. Moreover, we show that a harmonic trap can 
prevent the destructive effect of holes by confining them 
to the outside of the sample. Finally, via projected en- 
tangled pair states (PEPS) we complement our 
analysis with a simulation of the two-dimensional t — J 
model of hardcore bosons with antiferromagnetic inter- 
action. This setting is easier to investigate numerically 
than the corresponding 2D fermionic system and provides 
evidence that the physics studied in the one-dimensional 
case can be extrapolated to understand the conditions of 
an equivalent scheme in 2D. 

In the following, unless stated otherwise, we will focus 
on the one-dimensional case and consider a two compo- 
nent Fermi gas in an optical lattice. The physical set- 
ting consists of two adiabatic stages, depicted in Fig. [T] 
The first transition has already been realized experimen- 
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FIG. 1: The proposed adiabatic protocol. First, a BI in a lat- 
tice with depth Vi (a) is transformed to a product of singlets 
(b) by slowly switching on a second lattice with depth 
V2 and half the original wave length. Then, by lowering the 
barrier Vi —J- 0, the system turns into an AFM (c). 



tally |14| - |16| and can be straightforwardly described, so 



that we can focus on the second one and take Fig. [TJb) 
as the initial state for our theoretical study. In this 
situation, the system is governed by a one-band t — J 
model. This model emerges in the limit of strong inter- 
actions from the Hubbard model, which describes ultra- 



cold atoms in optical lattices [17| • We consider a bipartite 
t — J Hamiltonian with different couplings for even and 
odd links, H = Hp + Ho, where 

Hi = -h J2 (cI,,Cfc+i.,+H.c.) (1) 

rikrik+i ' 



SfcS 



£ = e, o. 



The superexchange interaction Ji and the tunneling pa- 
rameter ti are related through the on-site interaction 
U as Ji ^At'^/U. We fix the couplings on the even 
links, te — t and Je = J, and choose a linear ramp- 
ing of the superexchange interaction on the odd links, 
over total ramping time T so that Jo(t) = J ■ t/T and 
t„(r) = t ■ ^JVjT, for < T < T [13. In the following, 
we set J = 1 . A harmonic trap is included by adding a 
term Vt J2ki^ - ^o f^k to Eq. (1). 

The ramping time from Fig. \V[h) to [ijc) needs to be 
long enough such that the final state is close to the true 
AFM. The required time to ensure a certain degree of 
adiabaticity can be seen to scale as T oc 1/ , where 
A is the minimum gap between the ground and the first 
excited state during the evolution. A closer look at the 
relevant energy levels reveals that in this adiabatic tran- 
sition the gap decreases monotonically from J to the 
Heisenberg gap, which vanishes in the thermodynamic 
limit, and there is no phase transition occurring in be- 
tween 

The adiabaticity of the evolution and the fi- 
nal antifcrromagnctic order can be probed by 



two experimentally accessible observables. The 
first one, the squared staggered magnetization, 
MLs = W^E^m=ii-^y^"'i^i^rn), is the autifcr- 
romagnetic order parameter and can be determined by 
noise correlations [21[ . The second one is the double well 
singlet fraction, Po = jrY^kee il ^ i^kSk+i)) ■ Since 
the initial state. Fig. [Ijb), has a pure singlet in each 
double well, measuring Pq ^ 1 at the end indicates 
a change in the state. The generation and detection 
of singlet and triplet dimers in double well lattices 
has been recently reported [l^. Whereas the squared 
staggered magnetization is experimentally detected 
over the whole sample and captures information on 
long-range correlations, the singlet fraction can be 
determined in situ and hence allows us to probe parts 
of the sample. For the Heisenberg antiferromagnetic 
chain in the thermodynamic limit, these observables 



take on the values Afg^^^^g = and Pq.td ~ 0.693 
while in 2D M^^g ^pQ ss 0.0945 and Pq.td ~ 0.585 
Notice that for the finite ID systems considered in this 
work, Afjpjjjg does not vanish, but has a sizable value. 



comparable to the 2D thermodynamic limit [24 1. 

Absence of holes. Using the numerical simulation of 
the chain dynamics with MPS, we investigate the state 
at the end of the protocol for varying ramping time. To 
characterize the antiferromagnetic order independently 
of the system size, we define the relative quantities 



ML, 



iT)/ML,,AFM, PoiT) := Po{T)/Po,AFM, 



and espin(T) := £^spin(T')/i^spin,AFM, where the denomi- 
nator is the expectation value of the observable in the 
true AFM for a given lattice [23|. For the last quantity, 
-Espin is the expectation value of the spin term in the total 
Hamiltonian of Eq. (1). 

Figure O shows our results for an ideal case with no 
holes in the initial state, with all relative quantities 
converging to 1, as expected, in the limit of large T 
[Figs. [Ua),(b)]. The ramping time necessary to reach a 
certain relative magnetization grows with the system 
size. If we study, given T, which is the largest system for 
which a fixed value can be achieved [Fig. [UJa) inset] , 
we find N'^ cx T, which is consistent with the adiabaticity 
condition for a gap closing like A oc 1/7V [2^. 

For very long chains, the required T might not be ex- 
perimentally feasible. Remarkably, this does not exclude 
the preparation of antiferromagnetic order on large sys- 
tems. We may evaluate the magnetization over a sublat- 
tice in the center of the sample. If, given T and N, we ask 
for the largest sublattice size L for which the magneti- 
zation reaches a fixed value [Fig. Eljc)], we find a scaling 
T (X L?, as governed by an effective local gap, and not 
by the gap of the total system. In contrast to m? , the 
observables po and Cgpin do not depend on L [Fig. [UJd)] . 
This can be understood from the fact that po and Cgpin 
are determined by a two-site observable S^Sfc+i averaged 
over a sublattice of length L, whereas w? is a true L-site 
observable, and thus effectively probes the adiabaticity 
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FIG. 2: Absence of holes, (a), (b) m^, po, and Cspin as func- 
tions of the ramping time T, for chain length N = 22 (solid 
blue lines), 42 (dashed red lines), 62 (dash-dotted green lines), 
and 82 (dash double-dotted brown lines). The inset of (a) 
shows the squared size A'^'^ of the longest chain reaching a 
fixed = 0.85 at ramping time T, and reveals the scaling 
T oc iV^ (c), (d) Same quantities as above, evaluated on sub- 
lattices of length L = 22 (solid blue lines), 42 (dashed red 
lines), and 62 (dash-dotted green lines) for = 82. Now, the 
inset of (c) shows the squared size of the largest sublat- 
tice reaching = 0.85 at ramping time T, and reveals the 
scaling T oc . All results were obtained with MPS of bond 
dimension D = 60 and Trotter step St = 0.02 [H]. 



on the sublattice. These observations can be addressed 
analytically by means of spin wave theory (27j . The ex- 
perimental consequence is that with a large sample, high 
values of can be obtained in short ramping times T 
on small parts of the system, L cx a/T. 

Effect of holes. In a real experiment, the finite tem- 
perature causes the sample to be in a thermal mixture. 
As a consequence, localized holes will be present in the 
initial state, Fig. [Ub). Since the double wells arc de- 
coupled, the wave function of a hole will be an equal su- 
perposition of being in the left and in the right side of a 
single double well. Our simulation reveals that holes have 
a highly destructive effect on magnetic order. As seen in 
Fig. [3l^a), a few holes initially located on the boundary 
of the sample are enough to cause a dramatic reduction 
of the final staggered magnetization. 

We observe that the dynamics of holes can be qual- 
itatively well understood using a simplified picture, in 
which the spreading of an initially localized hole, prop- 
agating in an antiferromagnetic background, is modeled 
by a free particle. This picture is accurate in the limit 
i » J, when the spin term in Eq. (1) is negligible. We 
checked that it is also valid in the whole experimentally 
reasonable parameter regime by comparing the behavior 
of single holes to that of free particles with the same ini- 
tial wave function for the relevant range of times 
and tunneling values, 2 < t < 8, the hole spreads like a 
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FIG. 3: Effect of holes and harmonic trap, (a) rri^ as a 
function of the ramping time T, evaluated on a sublattice of 
length L = 42 for TV = 86, without holes (thin solid line), and 
with initially 2 holes at each boundary (thick solid line), and 
t = 2. The inset shows the largest sublattice size L reaching 
m? = 0.85 at ramping time T, for the two cases of the main 
plot, and the size of the hole-free region (dashed red line), 
(b) m? evaluated on a middle sublattice of length L = 82 for 
N = 102, with no holes (thin solid line), and with initially 
10 holes at each boundary and a harmonic trap of strength 
Vt = 0.004 (dash-dotted green line), 0.006 (dashed red line), 
and 0.02 (solid blue line), and t = 3. For Vt = 0.2 (not shown) 
the exact behavior of the ideal case is recovered. Again, the 
inset shows the largest sublattice size L reaching = 0.85 
at ramping time T, for the cases of the main plot. All re- 
sults were obtained with MPS of bond dimension D = 60 and 
Trotter step St = 0.02 [5^]. 



free particle with a maximal velocity v = 2t. In the case 
of small i = 2, a hole with higher initial energy causes a 
higher spin excitation, while for a large t = 8, the hole 
excites the spin background by AiJspin ~ 0.5 in the be- 
ginning, independent of its initial kinetic energy. This 
can be understood by assuming a simple classical Neel 
background, where the dclocalization of a hole, initially 
positioned at a boundary, breaks up exactly one anti- 
ferromagnetic bond. This assumption should become a 
good approximation for the hole dynamics in the regime 
t 3> J, where the timescale of the dclocalization is much 
faster than the timescale of the reacting spin background. 
In all cases, the squared staggered magnetization is re- 
duced substantially during propagation of the hole until 
it reaches a minimum after the hole has travelled once 
over the whole sample. We found that the magnetization 
reduction depends only weakly on the initial kinetic en- 
ergy, and on but it depends strongly on the number of 
holes in the sample. 

The simplified free particle picture allows us to inter- 
pret the results from Fig. [HJa). In particular, the strong 
magnetization drop around T « 8 indicates the arrival of 
holes at the middle sublattice. We can roughly predict 
this arrival time from the spreading of a free particle wave 
function, which, after time T, will have covered |t • T 
sites, taking the ramping of the lattice into account. We 
observe that the magnetization for a given sublattice L 
behaves like in Fig. [2][c) only for short ramping times, 
while the region is hole-free, until the holes reach the 
sublattice [Fig. [2Ia) , inset] . 
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Harmonic trap. The simplified picture described 
above points out that the negative effect of holes can be 
controlled by the presence of a trap. An external poten- 
tial changes sign for a hole and effectively turns into an 
inverse trap, capable of confining the holes to the outer 
parts of the chain. The trap strength should be chosen 
as large as possible without exceeding the on-site inter- 
action U, what would destroy the MI Fig. [TJc). 

The results of the dynamics within the harmonic trap 
are shown in Fig. |3|b). From energy considerations, a 
hole delocalizes at most by ±2i 2j|. As the trap strength 
is increased, the holes get more localized on the outside. 
As a consequence, the magnetization of the total sam- 
ple increases, and the behavior of the hole-free case is 
recovered. 

Two-dimensional case. For the 2D case, the adiabatic 
setting consists of an array of initially decoupled chains 
like Fig.[ljb), connected by a transverse lattice with the 
time-dependent couplings Jo and to- Different to the ID 
case, this system exhibits a phase transition in the ther- 
modynamic limit at Jo ~ 0.5 (20j . 

Although the numerical simulation of the 2D setting 
is much more demanding than the one for chains, we 
can relatively easily obtain results for hardcore bosons 
on lattices of moderate size. In the absence of holes, the 
bosonic and fermionic t ~ J models are equivalent, and 
our simulations serve to test our protocol on a 2D system. 
Similarities between both models under inclusion of holes 
are a subject of current research [2^, but our simulations 
can still provide a qualitative indication of the controlling 
effect of the trap. It is worth noticing that the AFM can 
also be realized with ultracold bosons [2^ . 

In the ideal case of no holes, we observe [Fig. lU^a)] 
that, whereas the energy converges quickly, the magneti- 
zation does not, and thus we cannot claim convergence to 
the true AFM within the numerically accessible ramping 
times studied here. Remarkably enough, we find that for 
a 10 X 10 lattice, a significant magnetization value m? is 
obtained in times of the same order of magnitude as for 
a chain of length 10, suggesting that the generation of 
antiferromagnetic order on much larger 2D lattices will 
be experimentally possible within reasonable timescales. 
Upon hole injection, a similarly dramatic magnetization 
reduction is observed [Fig. IH^b)] which can be controlled 
by the presence of a harmonic trap, as in the ID case. 

Discussion. We have proposed and analyzed an adi- 
abatic protocol, suitable to prepare an antiferromagnet- 
ically ordered state in an optical lattice, even from an 
initial state containing defects. The timescales for the fi- 
nite systems studied in this work lie well within the range 
of current experiments. Furthermore, we have observed 
that antiferromagnetic order can be produced in a sub- 
lattice in times governed only by its size. 

This scheme offers several advantages over other pro- 
posals. First, starting from a BI simplifies the prepara- 
tion of a sufficiently low entropy initial state. Addition- 




FIG. 4: Performance in 2D. (a) and egpin as functions 
of the ramping time T, for systems of size A'' = 4 x 4 (solid 
blue line), 6x6 (dashed red line), 8x8 (dash-dotted green 
line), and 10 x 10 (dash double-dotted brown line), (b) 
for A'^ = 8 X 8, with initially 4 holes at the boundary and 
no harmonic trap (dash-dotted green line), and with a trap 
of strength Vt = 0.25 (dashed red line), and 2.5 (solid blue 
line), and t = 2.5. The results were obtained with PEPS of 
bond dimension D = 4 (a) and D — 2 (b) and Trotter step 

St = 0.03 m. 



ally, the initial ground state in Fig. [Ijb) already features 
the final SU(2) symmetry of the AFM, so that the num- 
ber of excited states to which the evolution couples is 
minimum, as compared with an alternative proposal Q 
with only U(l) initial symmetry. Moreover, since hole 
doping is experimentally feasible, the same procedure 
can possibly be used to prepare the ground state with 
varying hole densities. This would open the door to the 
experimental exploration of open questions in condensed 
matter theory, ultimately the existence of d-wave super- 
conductivity in the t — J model. 
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Supplementary Material 

111 this supplementary material, we describe the numerical method used in our article, provide absolute values of 
experimental observables and discuss the numerical errors. 

Numerical method 

We employ Matrix Product States (MPS) in ID and Projected Entangled Pair States (PEPS) in 2D to compute 
ground states and time evolution [H-Q- For a chain of N quantum systems, each having physical dimension d, a MPS 
has the form 

l^MPs) - tr(A[l]«iA[2]^-...A[iV]^~)|siS2...Siv) 

SI .32 SN 

where the Si run from 1 to c? and the ^[1]"', for a fixed Si, are matrices of dimension D x D. The number of 
variational parameters in the ansatz is determined by the bond dimension, D, which also bounds the entanglement 
of the MPS. The MPS family constitutes a very good approximation to ground states of gapped local Hamiltoiiians 
in one dimension Q, and it has become a very successful tool for the study of quantum many-body systems. PEPS 
constitute the natural generalization of the MPS ansatz to larger dimensions [5|, |6|, and they are also known to be 
good approximations for ground and thermal states of gapped local Hamiltonians 18| . 

The basic algorithms for the study of many-body systems using these families of states are the variational search for 
the ground state and the simulation of time evolution. A detailed description of these algorithms, which are similar for 
ID and 2D systems, can be found in the review paper (see also references therein). The computational cost of MPS 
algorithms scales as 0{ND^), which allows the simulation of large chains and bond dimensions. PEPS algorithms 
scale as 0{ND^^), so that 2D simulations arc much more demanding, and only smaller values of the virtual bond 
dimension D can be considered. 

In the variational search for ground states, the numerical error comes from limiting the bond dimension to a 
certain maximum value. Comparing the results with those obtained after running the search with a larger D gives an 
estimation of the magnitude of this truncation error. Typically, the algorithm is run repeatedly with increasing bond 
dimension until convergence is achieved within the desired numerical precision. 

The simulation of time evolution of MPS (or PEPS) is based on a Suzuki- Trotter decomposition of the evolution 
operator This introduces another source of numerical error, in addition to the truncation of the bond dimension. 
The magnitude of this Trotter error can be controlled by decreasing the size of the time step, St, or using a higher 
order decomposition of the exponential. In particular, we use a second order Trotter decomposition in the case of one- 
dimensional simulations, while a first order decomposition is used in 2D. For each time step, the evolution operator 
for 6t is applied and an optimal MPS or PEPS approximation to the evolved state is found as described in [J]. The 
magnitude of the Trotter error is controlled by comparing results for various values of St, and the truncation error, 
as in the ground state search, is estimated from the comparison of results for different D. 

Absolute values of experimental observables 

The relative quantities shown in the main text result from normalizing the computed expectation values at the 
end of the evolution to the corresponding values in the true AFM ground state. For completeness, we show in this 
section the actual absolute values obtained for each observable, as well as the reference AFM values, and discuss the 
convergence of the numerical results. 

As explained in the article, we model our system by a bipartite t — J Hamiltonian, that, in ID, reads: 

H = -tc E (4,.c,+i.. +h.c.) + JeE(SfcSfc+i - ^^^)-to E (4^,c,+i,. +h.c.) + ,/oE(SfcS,+i - ^^^) ,(2) 

where the subscripts e and o denote even and odd sites and double occupancy of sites is forbidden, as implicitly 
assumed for the t — J model. The couplings on even links are constant, t^ = t and Jo = J, while the time-dependent 
odd couplings are increased from to their final values, t and J, during a total ramping time T, according to 
to{T) = t ■ yr/T and Jo{t) = J ■ t/T. We set J = 1. The two-dimensional system consists of several such chains, 
coupled in the transverse direction by to{T), Jo{t). 



7 



N 


MLgiD = 80) 


MLgiD = 100) 


PoiD = 80) 


Po{D = 100) 


E,pin{D = 80) 


E,pin{D = 100) 


22 


0.139654860817 


0.139654869116 


0.807121935054 


0.807121935023 


-0.434912539817 


-0.434912539817 


42 


0.086355374176 


0.086355479294 


0.775723907194 


0.775723905709 


-0.438751678570 


-0.438751678582 


62 


0.064239765116 


0.064240094937 


0.760657976140 


0.760657965104 


-0.440148664874 


-0.440148664972 


82 


0.051802232654 


0.051803243348 


0.751408113523 


0.751408046065 


-0.440871682329 


-0.440871682728 



TABLE I: Squared staggered magnetization M^^^^, mean singlet fraction per double well Po, and mean spin energy per site 
-Espin, for the AFM of total length A'^ obtained from ground state computation with MPS of bond dimension _D = 80 and 
D = 100. 



L 


MLs{D = 80) 


AfeLg(I?= 100) 


22 


0.141157102848 


0.141157444444 


42 


0.088248384680 


0.088248983897 


62 


0.065508480733 


0.065509337169 



TABLE IL Squared staggered magnetization A/ft^g for middle sublattices of length L — 22, 42, and 62, for an AFM of total 
length A'^ = 82 obtained from ground state computation with MPS of bond dimension _D = 80 and D — 100. 



Reference values in the AFM ground state 



We use the algorithms described above to compute numerically the true AFM ground state for various lattices, as 
shown in Tab. |T] (for chains of lengths iV = 22 — 82). The numerical convergence is checked by comparing the values 
obtained using bond dimensions 13 = 80 and 100. As can be seen from the values in the table, the maximum relative 
error is of the order 10~^ for the magnetization. Note that the energy is correct up to 10~^. We use these values as 
the reference for adiabaticity of the total chain of length N. 

As discussed in the main text, for long enough chains, we observe that antiferromagnctic order develops in a middle 
sublattice faster than on the total chain. We find that the timescales for observables measured on the sublattice are 
controlled by the range of the observable itself, as far as finite size effects can be ignored. Therefore, in order to quantify 
this observation, we need to compare the observables in the evolved sublattice with the corresponding AFM values for 
a sublattice of the same size in an infinite chain. The thermodynamic mean singlet fraction Po,td = lii(2) « 0.693 and 
mean spin energy E'spin.TD = 1/4 — ln(2) « —0.44315 are well known [l3|- The squared staggered magnetization on a 
finite sublattice, however, cannot be computed exactly, so that we approximate the reference value by the numerical 
estimation in a long chain (N = 82), shown in Tab.HIl Increasing the chain length, the reference values do not change 
significantly, as one can see in Tab. Illll for N = 142. 

In 2D, the corresponding values of the AFM were obtained with Quantum Monte Carlo by means of the ALPS 
code [ll|, and they are listed in Tab. II VI 



Adiabatically evolved states 



In the following, we present the absolute values of our observables for the state obtained at the end of the adiabatic 
ramping, in the same sequence as they appear in the main text. 



L 


MLsiD = 80) 


MLg(Z3 = 100) 


22 


0.141068047422 


0.141069521513 


42 


0.088395582256 


0.088397878145 


62 


0.065999995175 


0.066003286606 



TABLE III: Squared staggered magnetization M^mg for middle sublattices of length L = 22, 42, and 62, for an AFM of total 
length N = 142 obtained from ground state computation with MPS of bond dimension D = 80 and D = 100. 
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N 




-f/spin 


4 


X 


4 


0.233090090642115 


-0.574325441574560 


6 


X 


6 


0.164620512963000 


-0.603311944444444 


8 


X 


8 


0.135035306250000 


-0.618890781250000 


10 


X 


10 


0.121902426000000 


-0.625150439774003 



TABLE IV: Squared staggered magnetization M^t^„ and mean spin energy per site i?spin for the AFM of total size TV obtained 
from Quantum Monte Carlo using the ALPS code [lH . 



Absence of holes 

Fig. [5] and |6] show the computed expectation values at the end of the protocol, as a function of the total ramping 
time, for the case of no holes. Convergence of the numerical results is checked by comparing the results for bond 
dimension D = 40 and 60 and for Trotter steps 5t = 0.02 and 0.005. The largest relative error found (for the largest 
system and the longest ramping time) is of the order of 10~^, which ensures enough precision for our analysis. 




10 20 30 40 50 10 20 30 40 50 10 20 30 40 50 



T T T 

FIG. 5: A£ftag, Poi and -Bspin, as functions of the ramping time T, for chains of length N = 22 (solid blue), 42 (dashed red), 
62 (dash-dotted green), and 82 (dash double-dotted brown). The results were obtained with MPS of bond dimension D = 60 
and Trotter step St = 0.02 (lines), D = 40 and St = 0.02 (circles), L> = 40 and St = 0.005 (crosses), and L» = 60 and St = 0.005 
(squares) . 




10 20 30 40 50 10 20 30 40 50 10 20 30 40 50 
T T T 



FIG. 6: Same quantities as in Fig. [5] evaluated on sublattices of length L = 22 (solid blue), 42 (dashed red), and 62 (dash-dotted 
green), of a total lattice of length = 82. 



Effect of holes 



Injecting holes into the sample results in a substantial drop of the squared staggered magnetization, and an increase 
in the energy (Fig. [7]). The latter implies that the system gets excited and the numerical simulation via MPS becomes 
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more demanding. Now, the relative error after the longest ramping time, T = 30, for the case of 4 holes on A'' = 82 
sites shown in the main text, is of the order of 0.01 — 0.1, but it becomes significantly smaller for shorter times. This 
worst-case error does not affect our conclusions, since the main efi^ect we observe, the drop of the magnetization value 
upon hole arrival, is much larger (« 30%) and occurs already at much shorter times (T 20) , for which the numerical 
error is only of the order of 10"'^. The figure also shows that 2 holes on 42 sites are more dramatic than on 82 sites, 
and that the negative effect of holes increases with their number. 




10 20 30 10 20 30 10 20 30 



T T T 

FIG. 7: Ms^tag, -Poj and -Bspin, as functions of the ramping time T, evaluated on the middle L = 42 site sublattice, for 2 (solid) 
and 4 holes (dashed) on TV = 82 sites (thick blue) and = 42 sites (thin red), and the holes are initially located at the 
boundaries, and t = 2. The results correspond to D = 60 and Trotter step 5t — 0.02 (lines), _D = 40 and St = 0.02 (circles), 
D = 40 and St = 0.005 (crosses), and D = 60 and St = 0.005 (squares). 



Harmonic trap 



By including a harmonic trap Vt J2ki^ ~ koffik, the holes can be confined outside of the sample. We consider 10 
holes left and 10 holes right of a sample of size 82, and successively increase the trap strength, as shown in Fig. [51 
Consistent with the results in the previous section, our simulation is most demanding for the weakest trap, when holes 
can still enter the sample and excite the system. The largest relative error in that case is of the order of 0.01 — 0.1, 
but it decreases significantly if the trap strength is increased. Again, this worst-case error does not affect any of our 
conclusions. 




10 20 30 10 20 30 10 20 30 



T T T 

FIG. 8: Mg^tag) and -Bspin, as functions of the ramping time T, evaluated on the middle L = 82 site sublattice, for 10 holes 
initially on each boundary of 82 fermions, with a harmonic trap of strength Vt = 0.004 (dash-dotted green), Vt = 0.006 (dashed 
red), and Vt = 0.02 (solid blue), and t = 3. The inset in b) shows the occupation n of lattice site I after ramping time T = 30, 
and we find that the holes delocalize precisely ±2t at the boundaries of the trap. Again, the results correspond to D = 60 and 
Trotter step St = 0.02 (lines), L» = 40 and St = 0.02 (circles), L» = 40 and St = 0.005 (crosses), and D = 60 and 5t = 0.005 
(squares). 
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Two dimensional case 



In 2D, the time evolution is done with PEPS of bond dimension D ^ 2, D ~ 'i and D ~ A, and the largest relative 
error is of the order of 10^^, for the largest system in Fig.[9l Our results suggest that qualitative insight can already 
be gained from PEPS with D = 2. Therefore, Fig. [TOl shows the effect of holes and harmonic trap for Z3 = 2 without 
convergence check. Just as in ID [Fig. Elja)], the staggered magnetization decreases significantly with increasing 
number of holes, and the harmonic trap confines the holes on the outside. 




FIG. 9: A'fg^tag, Po, and -Bspin, as functions of the ramping time T, for A'^ = 4 x 4 (solid blue), 6x6 (dashed red), 8x8 
(dash-dotted green), and 10 x 10 (dash double-dotted brown). The results were obtained with PEPS of bond dimension D = 4 
and Trotter step 5t — 0.03 (lines), D — 2 (circles), and D = 3 (crosses). 




FIG. 10: (a) Mg^tag as a function of the ramping time T, for A'^ = 8 x 8 without holes (solid blue), with 1 hole (dashed red), 
and 2 holes (dash-dotted green), where the holes are initially localized at the boundary, and t — 2.5. (b) A'^ = 8 x 8 with 4 
holes initially distributed at the boundary, with no harmonic trap (dash-dotted green), and with a trap of strength Vt = 0.25 
(dashed red), and K = 2.5 (solid blue), and again t — 2.5. The results correspond to _D = 2 and Trotter step St — 0.03. 
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